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Abstract 

Let ffl, M and W be the free additive, free multiplicative, and boolean additive convo- 
lutions, respectively. For a probability measure ji on [0,c») with finite second moment, 
we find the scaling limit of (^^^^H-'v as N goes to infinity. The 7?,-transform of the 
limit distribution can be represented by the Lambert's W function. We also find similar 
limit theorem by replacing the free additive convolution with the boolean convolution. 

1 Introduction 

In free probability theory, some limit theorems are known as in classical probability theory. 

The most famous limit theorem is the free central limit theorem, which was found by 

Voiculescu. If {XjjjgN is a sequence of freely independent identically distributed (for short 

X\ "1- * * * A'jY 

free iid) random variables with finite second moment, then the normalized sum — 



'N 

converges to the Wigner's semi-circle law in distribution as N goes to infinity. In addition, 
we know the Poisson limit theorem, the stable limit theorem and so on, for details, see J9|, 
[5|], [3|]. Recently other new U,n. theorems „.h respect to the free convolutions fl. Q 
and [IJ] have been investigated. 

In this paper, we consider not only addition but also multiplication of free iid random 
variables. Especially, we introduce a new normalized sum of multiplications of freely in- 
dependent random variables, that is, for the double sequence of free iid random variables 
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having distribution on [0,oo) with finite second moment we consider a 
new normahzation Y/v 



where m\ is the mean of distribution ^. We shall see that the limit distribution depends only 
on the first and second moments. In its proof, we shall investigate the Taylor type expansion 
of the S'-transform. A similar limit theorem can be found under boolean independence. 

In Section 2, we shall gather the tools for free and boolean probabilty. Especially, we 
recall TZ, S and S-transforms and infinite divisibility. In Section 3, we shall prove the Taylor 
type expansions for S and S-transforms and some limit theorems. Finally, in Section 4, we 
discuss the limit distributions with focusing on infinite divisibility and moments. 

2 Preliminaries 

Let M+ be [0, +oo) and C+ be {z = x + iy G C; y > 0}. We fix notation that V and V+ mean 
the set of all Borel probability measures on M and respectively. We denote the free 
additive, free multiplicative and boolean additive convolutions by ffl, Kl and tt), respectively, 
see for details of convolutions, [3], [3] and [lit. Afterword, 5o stands the Dirac probability 
measure concentrated on 0. 

2.1 Analytic tools 

Here, we will gather the analytic tools for free and boolean probability and mention some 
their important facts. 

We denote the Cauchy transform of a probability measure ^ on R by 



and 



G^{z) = [ —^iJ.{dx), z e 
JrZ-x 



1 — xz 



denotes the moment generating function of p on R+. Then the Speicher's R and Voiculescu's 
7^~transforms of /x are defined as follow: 

R^{z) = zn^{z) = zG-\z) - 1, l/z e r„,;3, 

where G^^(z) is the inverse of G^{z) with respect to composition and r^^^^ = {z = x + iy £ 
C^;y > /3, \y\ > ax}. Note that we use both R and 7^~transforms for convenience. The S 
and S-transforms of p are defined by 

(z + l)^-i(z) 
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and 



respectively, where '^p^{z) is the inverse of ^'p(z) with respect to composition. Now, we 
summarize the relations between the transforms and convolutions, see for proofs 0] and [ll. 

Proposition 1. For fii ^V, ^2 ^ , Pi £ 'P+ oind p2 € V+, which are not 5o, there exist 
a > and /3 > such that 

For c > 0, the dilation operator Dc on is defined by 

for any Borel set B on R_|_, where — B = {x G M; —x G -B}. In the paper they showed 

c 

Sd,{p){z) = ^Sf,{z), 

and 

A probability measure p is ffl-infinitely divisible if for any n G N there exists /i„ G P such 
that 

A* = /^n ffl • • • ffl /^n ■ 

^ V ' 

n times 

We denote the class of all ffl-infinitely divisible distributions by I^. 

Remark 2. We can define other infinite divisibility replacing ffl by Kl or l±). But for boolean 
convolution, all probability measures are l±)~infinitely divisible. So we shall not discuss tt)- 
divisibility any longer. 

The next proposition characterizes the ffl-infinitely divisible laws. 

Proposition 3. The foUowings are equivalent: 
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(1) /.G/ffl. 

(2) TZ^ has an analytic extension defined on C~ with value U R. 

(3) There exist unique 6^ € M and finite measure such that 

n^{z)=b^+ [ -^i^^idt), ZGC-. 

JR >- — tZ 

The above expression is called ffl-Levy-Khintchine representation, or simply Levy- 
Khintchine representation. 

Example 4. The typical examples of ffl-infinitely divisible distribution are the Wigner's 
semi-circle law, the Dirac's delta distribution and the free Poisson distribution irt with 
parameter t >0 having density 

TTt{dx) = max(0, (1 - t))6o{dx) + ^V^t - {x - 1 - i)2l[(i_^)2,(i+^)2] (x)dx. (2.1) 

The Levy measure and of the semi-circle law are 6q and 0, and the free Poisson law 
TTt has bfj, = t and i^^ = t6i. We put vri by vr. 

The following functional equation of the R and S'-transforms can be found in, for in- 
stance, 



12l . Lemma 2]: 



Proposition 5. Assume that /x G P^. For some sufficiently small e > 0, we have a region 
Df, that includes {—it; < t < e} such that 

z = {zS^{z)) , (2.2) 

for z £ Ds- 

3 New limit theorems 

In this section, we prove a new limit theorem related to both free additive and multiplicative 
convolutions. We also discuss a similar result with replacing ffl by tt). It was proved in (lo| 
by Mlotkwski that for the free Poisson law tt, we have 

i^n((vr^("-^))'")"^-0 indist., 



pP 

where the p^^ moment of uq is given by — . We find that theorem of this type holds for any 

p\ 

probability distributions under finite second moment condition. 
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3.1 Taylor type expansion 

We prove the Taylor type expansion for the 5-transform under the finite moment condition. 
Concerning with the 7^-transform, it was proved by Benaych-Georges in For each region 
A in C, we denote z — )• with z £ Ahy z — > 0. 

Lemma 6. Let p G Vj^ has the moment of order p, that is, for k = 0,1,2, ... ,p, 

mk{p) := / x'^ p{dx) < oo. 



Then its moment generating function ^ p{z) has a Taylor expansion 

p-i 



^p{z) = ^Mp)^'' + 0{zn, z 0. 



mk[p)z-\-U(z''), z 

k=l 

Proof. We have 



k=l 

Note that for given x G M+, the inequahty 

xPzP 



C x^z^ 

pi^) -Y,^k{p)z^ = / -. p(d, 



x). 



< \xz\ 



1 — xz 

holds for z with sufficiently small l^j. By the Lebesgue dominated convergence theorem, 

p-i 



lim ^ 1 ^p{z) -Ymk{p)z^ J = mp(p). 



.^^0- \ fc=l 

□ 

Lemma 7. Let p G "P^ has the moment of order p and p ^ 5q. Then we have the followings: 

(1) ^ p{z) is univalent in iC"*". 

(2) The inverse function : ^'p(iC~^) — )■ iC~^ of'^p admits the Taylor expansion 

V 



%\z) = Y,bk{p)z^ + o{zP), z'^'O. 



k=l 



In particular, bi = -— and 62 ^ ^ 
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(3) Dp := ^p(iC'*") is a region contained in the circle with diameter ( p({0}) — 1, O). In 
addition, ^p(iC+) nM = (p({0}) - 1, O), 



\im^-^(t) = 



lim ^^^(t)=oo. 

4P{{0})-1 P 



and 



Proof. (1) and (3) are proved in Bercovici and Voiculescu [5|, Proposition 6.2]. 

(2) Step 1 We shall first prove for p = 1. Take any continuous path {z{t)}te{o,i] in Dp 
such that limt^Q z{t) = 0. By (1), we can choose a unique continuous path {w(t)}(g[o^i] on 
t G (0, 1] such that lim4ow(t) = and ^p{u{t)) = z{t). 

*p'(^(i)) uj{t) 1 1 
lim — = iini = lim = . 

40 z{t) 40 ^p{uj{t)) 40 ^p{uj{t))/uo{t) mi 

By arbitarity of the paths, it follows that 

^-p\z) = hip)z + oiz), z'^'O. 

Step 2 Assume that the statement of (2) holds for I < p. From Lemma [6] we have 

I 



%{z) = miip) + J2i^' + l)mfc+i(p)z'= + o{z'), z 0. 



k=l 

By the assumption, we have 



Uk{IJ)Z ^ U{Z 

k=l 

Thus, 

I 

mk+i[p)z''' + o[z') z 

k=l 



o ^-\z) = + J2 f^k+i{p)z' + o{z') z 0, 

k=l 

where fh2{p) = 2m2{p){bi{p))'^ . Since the derivation of ^^""^(z) is given by 



p 

1 

we obtain 

I 

mi{p) 



fc=l 
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which imphes that 



k=l 

Comparing the first two coefficients, it follows that 

h{p) = — ^ 
mi{p) 

and _ 

OA / ^ "^2(/>) _m2{p) 
2b2{p) = ^ = -2- -—3, 

hence, 

As a consequence, we can derive the following theorem: 



□ 



Theorem 8. Let p G V+ has the moment of order p and p ^ Sq. We denote the variance 
of p by Var(p). 

(1) The S -transform of p admits the Taylor expansion of the form 

Sp{z) = ^Sk{p)z''+' + o{zn, z'^'O. 

k=0 

Futhermore, So{p) = and si{p) = - = 

(2) The Yj-transform of p admits the Taylor expansion of the form 

= J2'^k{p)z''+' + o{zn, z 0. 

fc=0 

Futhermore, ao{p) = l/mi(p) and ai{p) = ^ - = 

3.2 Limit theorems 

Here we state the main theorem. 

Theorem 9. We assume that p G V+ has the second moment and put a — —YsiM 



(1) There exist sq > and si < such that the S-transform of p is given by 

Sp{z) = So + siz + o(z), as z — f 0, 



7 



and there exists a probability measure t^Q, G V+ such that 

In addition, the S-transform of the limit distribution {)a is 

S^ciz) = exp(-az) . 
(2) There exist ctq > and ai < such that the T,-transform of p is given by 

T,p{z) = (To + aiz + o{z), as z 0, 
and there exists a probability TTiGusuTc Set G T^-j- sucfi that 

indist.. 

In addition, the T, -transform of the limit distribution 5a is 

^So.{z) = exp(-az) . 

Proof. By Propositon [H we obtain 

1 / 5l2^\ ^ 

(( mu\^^\ (2;) = — So H + 0(1) (for sufficiently large n) 

^sY^/n[\p j J Sq V n / 

n— >oo / Si 

— > exp — z 
\so 

From 01 and Theorem 6.13 (i)], there exists a free multiplicative infinitely divisible 

measure t)a such that S„ (z) = exp ( —z ) . The proof for (2) is the same as for the free 

\so J 

additive case. □ 

We can exchange the order of free multiplicative and free additive (or boolean additive) 
convolutions. The difference is in the scaling speed. 

Corollary 10. Under the same setting as in Theorem\^ we have 
as 77- —7- 00. 

Proof. As we have done in the proof of Theorem El it can be proved by using Proposition 
[Hand Theorem [8j □ 
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4 Lambert W function and infinite divisibility of the limit 
distribution 



4.1 On the limit distribution of free case 

When we calculate the i?-transform or the moment generating function, the Lambert's 
VF-function plays an important role. Let Wq{z) be the principal branch of the Lambert 
-function, that is, it satisfies the functional equation 

z = Wo{z)eMWo{z)). 

For more details of the Lambert W function, see, for instance, Q]- 

Theorem 11. (1) The TZ and R-transforms of probability measure rjo are given as fol- 
lows: 

-Woi-az) 



az 



1 



Rr,Az) = Woi-az). 



a 



(2) t)a is both ^-infinitely divisible and M-infinitely divisible. 



(3) The free cumulant sequence of \}a is 



[an 



in— 1 



ni 



nSN 



(4) The Levy measure v^^ of X}a is given by 

^^cids) = —sf~\a/s)lio^ac]{s)ds, 



aiT 



where f{u) = ucscuexp{—ucotu). In case of a = 1, for the shape of the density of 
u^^, see the graph below. 



0.5 



0.4 



0.3 - 



0.2 



0.1 - 



0.5 



1.0 



1.5 



2.0 



2.5 



(5) It holds the following formulas: 
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The proof of this theorem is helped by the following well-known integral representation 
of the Lamber'ts 14^-function: 

Proposition 12. For any z G C"*", we have an integral representation: 

Wojz) _ J_ r jl-ucotuf + u'^ 
z 71" Jo ^ + ticsc«exp(— ncot n) 

Proof of Theorem 11. (1) The Kl-infinitely divisibility is trivial from the form of the S- 
transform and the facts in 0]. By Propositon O we have 

R,^{ze-n = z. 

Then the i?-transform is given by using the Lambert's W^-function as follows: 

Rr,^^) = --Wo{-az), 
a 

and hence, 

Wo{-az) 



-az 



(2) By the property of Wq{z), TZ^^ has an analytic extension defined on with value 
U M, which means that t)a is ffl-infinitely divisible. 

(3) The Taylor type expansion of —Wq{—z) at the origin is obtained from Equation (3.1) 
of 0, pp. 339]. 

(4) We put g{u) = (1 -ucot{u))'^ + u^. Noting that 

= -J^, (4.1) 



we obtain 



' ' ttJo -az + f{u) 



1 r 9iu)/fiu) rrvzf). 



vr Jo ^ — ctz/f{u) a-K Jq 1 — sz 
where we have changed the variables as s = a/ f{u). 



air Jq 1 — sz 

= — r (t^ + i) r'Hs)ds 

a-K Jq \l-sz J 

1 1 /■"'^ 7 

= - + — sf-\a/s)ds. 

a an Jq 1 — sz 

sf~^(a/s) 

Therefore we obtain the Levy measure v„ ids) = ds of X)a- 

air 

(5) It is direct consequence of Proposition [TJ □ 
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Remark 13. Here we consider the limit distribution with paramete a = 1. For example, 
if p is the free Poisson distribution with parameter 1, this is the case. Simply we write t) 
instead of t^i. There exists a probability measure p such that 

7^p(z) = ^iM^. (4.2) 

Indeed, if we consider the shifted free cumulant sequence {kn(p)}neN = \ \ , 

which is a sequence of coefficients of Taylor expansion of i?p at 0, then it becomes a moment 
sequence of a probability measure. This means that the measure /) is a free compound 
Poisson distribution with a compound measure a, the moments of which are mn{(T) = 
■ From ()4.2p . we have 

zn„izMp{z)) = zMp{z). (4.3) 
By putting P{z) = zMp{z) and using the Lagrange inversion formula, (|4.3p implies that 
n*^ coefficient of {Piz)} = — x ((n - I)''* coefficient of 7^„(z)) . 



Hence we obtain the moments of p as mn{p) 



4.2 On the limit distribution in boolean case 

— - \ , the positivity 
n! J n>o 

of which is ensured by 10|. Then its moment generating function M^(z) can be given by 

M,{z) = y \ = - — — (1) 



where the function r]{z) is defined by 

T^^z) = -Wo{-z), zgC\[-^,oo). 



Remark 14. The following useful facts on the function rj can be found in [8|, Sect. 2]: The 
map 

exp {O cot 0) 



sin( 



IS 



a bijection of (0, vr) onto (0, e), and if we define the functions rj^ ,ri : [ — , oo) — )■ C by 



T]'^ exp ( - 61 cot 6*) = 6* cot 61 ± i 0, < 6^ < vr, 

V sin J 

then 

7y^(x) = lim ?7(x + iy), xg[-^,oo). 
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From (1), the Cauchy transform of the measure s is given by 

1 1 
V 



Ga(C) = - y-TT' forCGC\[0,e]. 



Now we apply the Stieltjes inversion formula to obtain the density function (ps{t) of the 
measure s, that is, for t G [0, e], 

ips{t) = limIm('Gs(t + ie] 

= lim Im 



TT eio y t + ie l — - 



t+ie 



lim Im i 



IT eiO I t2 _^e2 Y t-ie 

— Im I — 



where the function r]~ is defined as in Remark above. Here we change the variables 

1 



t sin I 



■ exp ( — 9 cot 9) , 



then it follows that 



1 T / 1 1 

^s[t) = Im 



vr \t I - {9coi9-i9) 
11 9 



^ t {l-9cot9) +9^ 
1/9 

■ exp { — 9 cot 9^ 



Vsin^^ ' ' ' ) \^(l-0cot 0)^ + 02^ 

1 9'^eyi^{-9cot9) 
^ sm9{{l-9coi9f + 9^ 

Thus we obtain the following proposition: 

Proposition 15. The probability density function ip^ of the measure s can be given by the 
implicit (parametric) form as 

sin V , X \ 1 exp ( — u cot v) 



V9g ( exp [y cot v) = r-, < u < vr 



^ sin v({\ — V cot v)^ + v"^ 



12 



Remark 16. (1) The shape of the density function of ip^ is as the graph below, especially 
non-unimodal. 



2.5 



2.0 




0.0 L 



0.0 



0.5 



1.0 



1.5 



2.0 



2.5 



(2) The function {\ — v cot v)'^ + v'^ also appears in the integral representation of VFo(2;) / z 



as we mentioned Propositonll2t 



[l — V cot v) + 



exp ( — cot v) 



dt). 



Thus using f{v) = f cscf exp(— iJcotu) and (|4.1|) again, the parametric form of the 
density function can be rewritten as 

v2 



( fiv 



1 {f{v)y 



f{v)J vr f>{v) 

Remark 17. We also know that there exists the similar moment sequence. In the paper 
by Dykema and Haagerup 0|, they find a limit distribution of DT-operator DT(0, Sq). The 
moment of their one is nin 



(n+l)! 
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